Application des modeles de Markov cache; 



1. Introduction : 

Apres avoir donne l'aspect theorique des methodes de Markov caches dans le chapitre precedent 
alors le moment est venu de les appliquer a un probleme reel, la poursuite d'une cible radar dans une 
zone prealablement determinee. 

2. Principe d'utilisation d'un HMM : 

II s'agit de determiner dans un modele prealablement appris la sequence d'etats correspondant a une 
observation donnee. On peut utiliser l'algorithme de [Baum Welch] de type "forward-backward" qui 
calcule la probabilite optimale d'avoir observe une observation O en sommant les probabilites 
associees a tous les chemins possibles permettant d'obtenir O dans le HMM. 

Mais l'algorithme le plus utilise est celui de Viterbi, il ne prend en compte que le meilleur chemin, 
c'est-a-dire celui qui maximise la probabilite d'obtenir O, l'algorithme de Viterbi est une solution 
recursive sous-optimale au probleme de l'estimation de la sequence d'etats dans un HMM, cet 
algorithme est egalement utilise dans le cadre plus general de la recherche d'un chemin sous- 
optimal dans un treillis. [12] 

3. Ressources : 

3.1 Materiel les : 

Pour elaborer notre travail nous avons eu besoin d'un Intel Pentium3 avec la configuration 
suivante : 

Intel Pentium 3 +1.00 GHz 
RAM 128 Mo 

3.2 Environnement : 

Systeme Sexploitation : Windows XP. 

Langage de programmation: Langage de programmation : MATLAB 6.5 
Logiciel qui permet de maniere interactive : 

> de faire des calculs matriciel; 

> d'analyser les donnees; 

> de visualiser les resultats en 2D et 3D. 

4. Probleme pose : 

Le probleme pose est le suivant: [10] 

Nous avons une zone bien determinees de 6x6 etats et 9 capteurs repartis de la maniere suivante 
Voir Figure 4.1. Ou chaque etat correspond a une position de la cible de coordonnees (m, n), 
avecl < m < 6, 1 < n < 6. 
Les 9 capteurs sont repartis dans la meme zone. 
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Chaque capteur est en mesure de detecter la cible en mouvement dans la position la plus proche. 




La legende: 

Represente un capteur dans notre cas c'est un Radar. 
: Represente les etats de 1 jusqu'a 36 respectivement de la gauche vers la droite. 




di : Represente la distance entre les etats (les positions de la cible). 
d2 : Represente la distance entre le capteur et la position de la cible. 

Dans notre travail on s'interesse a faire une experience qui conserve tous les capteurs ainsi que 
toutes les positions. 

Pour cela on peut traiter trois cas: 
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4.1 Le premier cas 




Figure 4.2 : Zone de travail des capteurs pour le premier cas 



Capteur 1 
Capteur 2 
Capteur 3 
Capteur 4 
Capteur 5 
Capteur 6 
Capteur 7 
Capteur 8 
Capteur 9 



il capte la position 1. 

il capte les positions 3 et 4. 

il capte la position 6. 

il capte les positions 13 et 19. 

il capte les positions 15 et 16 et 21 et 22. 

il capte les positions 18 et 24. 

il capte la position 3 1 . 

il capte les positions 33 et 34. 

il capte la position 36. 



On remarque que 16 positions sur 36 sont surveillees. 
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Figure 4.3 : Zone de travail des capteurs pour le deuxieme cas 



Capteur 1 
Capteur 2 
Capteur 3 
Capteur 4 
Capteur 5 
Capteur 6 
Capteur 7 
Capteur 8 
Capteur 9 



il capte les positions 1, 2 et 7. 

il capte les positions 2, 3, 4, 5, 9 et 10. 

il capte les positions 5, 6 et 12. 

il capte les positions 7, 13, 14, 19, 20 et 25. 

il capte les positions 9, 10, 14, 15, 16, 17, 20, 21, 22, 23, 27 et 28. 

il capte les positions 12, 17, 18, 23, 24 et 30. 

il capte les positions 25, 31 et 32. 

il capte les positions 27, 28, 32, 33, 34 et 35. 

il capte les positions 30, 35 et 36. 



On remarque que 32 positions sur 36 sont surveillees. 
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4.3 Le troisieme cas : 

Toutes les positions sont surveillees par les 9 capteurs. 

5. Presentation du programme : 

5.1 Constitution du corpus : 

La base de donnees est formee d'une matrice de transition (36x36) et d'une matrice d'observation 
(9x36). 

9 : represente le nombre de capteurs 

36 : represente le nombre des etats (i, j) position dans lesquelles se troue la cible. 
La longueur d'observation est fixee a 50 

5.2 Les fichiers MATLAB : 



Voici la liste du programme et des sous programmes qui seront executes sous MATLAB 6.5 sur la 
plateforme PC Windows: 

> HMM-application.m. 

> Matrice-transition.m. 

> Matrice-observation.m. 

> Gene-etat-vraie.m. 

> Viterbi.m. 

> Cumuprod.m 

> Calcul-erreur.m. 

Et maintenant nous presentons la methodologie de programme principale. 
5.3 Programme principal: HMM-application.m 

Ce programme est base sur les procedures suivantes: 

> Localisation et specification du nombre des etats. 

> Specifier la longueur de la sequence du temps pour produire l'observation. 

> Appel de la fonction Matrice-transition.m pour generer la matrice de probabilite de transition. 

> Appel de la fonction Matrice-observation.m pour determiner observabilite de chaque capteur. 

> Performer 20 iterations. Pour chaque iteration, on fait les etapes suivantes: 

• Appel de Gene-etat-vraie.m pour generer l'etat actuel de la cible. 

• Generer l'observation binaire associee avec chaque observateur, basee sur 
l'observabilite des capteurs et l'etat actuel de la cible. 

• Appel de la fonction Viterbi.m pour trouver l'etat estime. 

• Appel de la fonction Calcul-erreur.m pour calculer et enregistrer les erreurs de 
distance entre les etats actuels et les etats estimes. 

• Calcul de l'erreur moyenne entre les positions sur le temps pour l'iteration courante. 

> Calcul des moyennes est erreurs pour 20 iterations et les representees dans un graphe. 
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5.4 Apprentissage: 

Le but de l'apprentissage est de determiner les parametres qui constituent les modeles tels que la 
matrice de transition Ay et la matrice d'observation By. 

> La matrice de transition est determinee en fonction de la distance qui separe les etats qui est 
note di. 

> La matrice d'observation est determinee elle aussi en fonction de la distance qui separe le 
capteur et la position de la cible qui est notee 62- 

Cas 01 : Le capteur capte la position la plus proche seulement voir Figure 4.2. 

Cas 02 : Le capteur capte les deux positions les plus proches voir Figure 4.3. 

Cas 03 : Le capteur peut capter jusqu'a trois positions. 

6. Resultats, courbes et commentaires : 

Cas 01: pour le premier cas on a eu les resultats suivants: 
Etat actuel: 



13 13 20 32 (26) 27 35 34 34 34 28 28 28 27 33 33 26 20 (31) 25 25 
25 27 27 34 33 (33) (2$) 34 29 (35) 36 29 34 28 27 32 26 20 22 17 22 



9 16 23 (29) 24 24 22 22 
Etat estime: 

25 25 25 25 (26) 32 32 32 32 32 27 27 27 32 32 32 31 31 @) 31 31 
31 31 31 31 31 (33)(28) 35 35 (35) 35 35 35 35 35 36 36 36 36 29 29 



29 29 29 (29) 29 29 29 29 

(^) : Les positions entourees d'un cercle rouge ont ete identifie (etat actuelle = l'etat estime). 
La courbe des trajectoires actuel et estime est representee dans la figure 4.4. 
La courbe d'erreur quadratique moyenne est representee dans la figure 4.5. 
Remarque : 

La representation des courbes d'erreur quadratique moyenne dans les trois cas est basee sur le 
calcul de l'erreur moyenne sur 20 trajectoires distinctes. 
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Figure 4.5 : Courbe d'erreur pour le premier cas 
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Pour le premier cas (casl), on constate que l'erreur moyenne entre l'etat desiree et l'etat estimee est 
assez importante par ce que chaque capteur est capable de capte la position la plus proche 
seulement. 

Cas 2 : 

Dans ce cas, nous avons ajuste les parametres du modele qui sont en fonction des distances di et d2. 
On garde la meme matrice de transition du cas 1, et nous avons change la matrice d' observation qui 
est en fonction de la distance d2. C'est-a-dire on a varie d2de telle sorte que le capteur capte les deux 
positions les plus proches ce qui donne 32 position surveiller du totale 36 (voir la figure 4.3). 



Etat actuel: 

@ 3 10 10 22 33 28 28 27 @ 30 30 22 20 22 27 35 35 28 @ @ 
35 35 29 23 24 29 30 30 (36) 35 (36) ^6)(36) (36) 35 (35) (35) 24 24 (23) 17 



22 33 20 20 9 3 4 4 
Etat estime: 



0 4 4 4 4 4 16 28 28 @ 34 34 34 34 34 34 34 34 34 @ (34 



36 36 36 36 36 36 36 36 ^6j 36 (3J>) (3J>) (36) (36) 36 (35) (35) 23 23 (23) 23 
28 22 22 22 22 22 10 10 

La courbe des trajectoires actuel et estime est representee dans la figure 4.6. 
La courbe d'erreur quadratique moyenne est representee dans la figure 4.7. 



Cas 2 




Figure 4.6 : Courbe des trajectoires pour le deuxieme cas 
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Figure 4.7 : Courbe d'erreur pour le deuxieme cas 



Comme remarque nous avons constate une legere amelioration ou l'erreur moyenne entre l'etat 
desire et l'etat estime est importante. 
Cas 3 

Etat actuel: 

©©©©® 0 ©0® © ©© 1 1 2 2 2 2 ©0 2 

2 2 8 0 15 15 15 14 14 14 14 14 14 @ 19 19 @ @ @ @ @) 
0 19 @ 19 @ @ 7 8 

Etat estime: 

®®®®©®©©® © © © 2 2 3 3 3 3 ©©3 

3 3 9 0 9 9 9 8 13 13 13 13 13 @ 13 13 @ @ @ @ @ 
(13) 13 (13) 13 (l3) (13) 14 14 



La courbe des trajectoires actuel et estime est representee dans la figure 4.8. 
La courbe d'erreur quadratique moyenne est representee dans la figure 4.9. 
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Dans ce cas nous avons reajuste les parametres du modele. On a change la matrice de transition du 
cas 1 et cas 2 qui est en fonction de la distance di, et nous avons garde la matrice d'observation du 
cas 2 qui est en fonction de la distance d2. C'est-a-dire on a varie di de telle sorte que le capteur 
capte les trois positions les plus proche ce qui donne 36 positions surveiller et on obtenue des 
resultats satisfaisants (on constat une nette amelioration). 



La figure suivante represente l'erreur des trois cas en meme temps. 
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Figure 4.10 : Courbe d'erreurs des trois cas 



La figure 4.10 montre que le troisieme cas est le meilleur par rapport aux deux autre cas. 

Done on peut conclure que les meilleurs resultats peuvent etre obtenus par la diminution de la 
distance di. 



7. Conclusion 



Dans ce chapitre nous avons simule une application « La poursuite d'une cible RADAR dans une 
zone prealablement determinee » par la technique des modeles de Markov caches (HMM). 
La variation de la distance entre les etats d'une part et entre les capteurs et les etats d'autres part ; 
par l'utilisation de l'algorithme de Viterbi qui nous a permis des resultats satisfaisants par cette 
technique (figure 4.10). 

On peut dire que e'est un resultat encourageant pour ce genre de probleme, et que la technique 
(HMM) restera toujours valable dans plusieurs domaines d' application malgre leur ancienne 
apparition. 
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